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Abstract 

H. Yamamoto [1] has proposed employing B mesons produced in conjunc- 
tion with a single charged pion at an T resonance for studies of CP violation 
in the neutral B meson system at a symmetric e + -e~ collider. The sign of 
the charged pion would tag the neutral B meson. We estimate this branch- 
ing ratio, employing the heavy meson chiral effective field theory. We find a 
negligible branching ratio to BB^n^ at the T(5S) and a branching ratio of 
only a few percent at the T(6S). However, if nonresonant studies of neutral B 
mesons should prove feasible, Yamamoto's proposal could be a good method 
for tagging neutral S's for the study of CP violation at a symmetric collider. 

We also explore the possibility of studying B s at the T(5S). The rate is 
low but depends sensitively on the precise value of the mass of the B s . The 
background we compute is comparable to the rate at the largest allowed value 
of the B s mass. 

Finally, we discuss the extraction of the axial pion coupling to B mesons 
from measurement of the BBtt branching fraction in a restricted region of 
phase space, where chiral perturbation theory should work well. 
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1. Introduction 

It is well known that CP violation studies are difficult at a symmetric collider. 
In e + e~ — > B°B°, the heavy mesons are produced in a C-odd state, so that the time 
integrated asymmetry vanishes unless the time ordering of the signal and tag can be 
measured. An alternative tagging method for which the time integrated asymmetry does 
not vanish has been proposed by H. Yamamoto [1]. His suggestion is to study B mesons 
produced in conjunction with a single charged pion at an T resonance, so that the sign 
of the charged pion tags the single neutral B meson as B° or B°. Neutral B* mesons 
can also be used, since they decay immediately via photon emission into pseudoscalars, 
before weak mixing or decays have occurred. This method should provide a simple and 
efficient tag of the neutral B meson: In addition to the decay products of the neutral B, 
one only needs to detect the additional soft pion. The charged B is then tagged by the 
invariant mass of the missing four momentum, so it need not be reconstructed. Unlike 
many conventional proposals for the study of CP violation in the neutral B system, 
essentially all the events are tagged. 

The utility of this method depends on the event rate. We calculate this rate, 
employing the heavy meson chiral effective theory, in order to evaluate the potential of 
Yamamoto's proposal. Although it is difficult to do a reliable calculation for all relevant 
kinematic regions, we can nevertheless do a calculation which incorporates propagator 
enhancements, phase space, and derivative couplings. We find that the branching ratio 
should be negligible at the T(5S) and only a few percent at the T(6S). However, 
Yamamoto's proposal could prove a competitive method for tagging neutral S's for the 
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study of CP violation at a symmetric collider at slightly higher center of mass energy, 
about 12 GeV. 

Our calculation is also useful because the BBir mode is a potential background to 
B s identification (as discussed in section 8). We find the BBix background could be 
comparable to the rate for B s production at the largest experimentally allowed value of 
the B H mass, but is most likely small compared to the rate if the B s mass is near the 
central value of the experimentally allowed range. 

Finally, we discuss the extraction of the axial coupling constant of the B mesons 
from a measurement of the BBtx branching fraction in a restricted region of phase space, 
where the heavy meson chiral lagrangian should apply. 

We begin in section 2 by describing the heavy-meson effective theory, in order to 
review the assumptions and establish notation. In section 3, we describe how heavy 
quark ideas apply to the process e + e~ — > BB ± n Zf 1 and in section 4 we construct 
the effective lagrangian. We then calculate the cross section for e + e~ — > BB ± ir Zf in 
section 5. We discuss the implications of our calculation for CP violation studies and 
B s identification in sections 6 and 7. In section 8, we discuss the regime of validity of 
the results, and the extraction of g. Conclusions follow in the final section. 

2. Review of the Heavy Meson Theory 

In this section, we review the treatment of heavy meson fields in the heavy meson 
chiral effective theory. It is convenient to describe the B and B* mesons which contain 
a bottom quark of velocity v and a light antiquark, q by a Dirac tensor field of the form 
[2]: 

B(v) = 1 -±l(-b l5 + Eb s ?). (2.1) 
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Here b and b £ are the destruction operators for B and B* mesons. The field of the 
B mesons, which we denote as B(i>), can be obtained from B(t>) by using the charge- 
conjugation properties of the B mesons. One obtains [3] 

B» =C(CB(v)C- 1 ) T C T 

,1-4 ( 2 - 2 ) 

where C is the charge conjugation operator and C = i'j 2 ^ . 

We also require fields which create B mesons, obtained from the destruction fields 
by Dirac conjugation: 

B» = 7°B(v)t 7 

(2.3) 

B(v) = 7°B»V. 

Under a heavy-quark symmetry transformation the mesons transform as follows: 

B(v)^S v B(v), 

(2.4) 

B( V )->B( V )SJ, 

where is an element of the spin-1/2 representation of the little group of a particle of 
velocity v. 

We implement the chiral symmetry in the usual way [4] by introducing the non- 
linear field £ = e^ 7 "* 'I** . Under an SU(3)l x SU(3)r chiral transformation, £ — > 
L£f/t = U^B). This defines the SU(3) matrix U as a nonlinear function of L, R and 
7r(x). From £ one can construct an axial vector field, 

A » = Ltfd»z - e^t) = -^//, + • • • , (2.5) 

which transforms under the chiral symmetry as — > UA^U^ . One can also construct 
a vector field from £, 

F m = 1 (e t^ + ^t) = tt^/2/2 + . . . , (2.6) 



which functions as a connection term in a covariant derivative: under a chiral transfor- 
mation, + V* = -> UD^W. 

According to the usual convention [5], the field B a , which destroys a B meson 
containing a light antiquark of flavor a, transforms under the chiral symmetry as 



Similarly, the field B a which destroys a B meson containing an light quark of flavor a 
transforms as 



3. Applicability of the Chiral Heavy Meson Effective Theory to 
e+e - -> BB^* 

In this section we discuss the energy regimes in which B meson production can 
be reliably described using the chiral heavy meson effective theory. We argue that 
although such a calculation is not reliable to better than an order of magnitude within 
the resonance region, heavy quark methods should apply at higher center of mass energy. 
Because the calculation also requires that the pion energy be sufficiently low for chiral 
perturbation theory to work well, the calculation with the heavy meson chiral lagrangian 
is trustworthy only over a restricted region of phase space. 

In order to become familiar with the energy scales involved in BB and BBix pro- 
duction, we begin by listing the masses of the B mesons and the T resonances and 
also the total final-state kinetic energy, or Q value, for T decay into BB and BB^n^: 
The charged pseudoscalar B mesons have mass Mb = 5278.6 ± 2.0 MeV (all values for 



B a (v) ^ B a (v)Ul b . 



(2.7) 



B» ^U ab B b (v). 



(2.8) 
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Table 1. Masses and Q Values 



Resonance 


v^(MeV) 


Qbb (MeV) 


Qbb±,t (MeV) 


T(4S) 


10580 


23 




T(5S) 


10865 


308 


168 


T(6S) 


11020 


463 


323 


Off Resonance 


12000 


1443 


1303 



particle masses in this section are taken from ref. [6]); the neutral pseudoscalar mass is 

essentially the same. The vector B* meson mass is Mb* = 5324.6 ±2.1 MeV. In Table 1 

we list the central values for the masses of the T resonances and the Q values for decay 

into two heavy pseudoscalars with and without a charged pion: 

Q BB = M r -2M B , 

(3.1) 

Q BB ±^ =M r - 2Mb ~ m w ± . 
For final states including heavy vector mesons, the Q values are reduced by one or two 

times the heavy meson hyperfine splitting, Mb* — Mb = 46.0 ± 0.6 MeV. We also list 

the Q values for a center of mass energy of 12 GeV, which is above the resonance regime. 

We see that pions produced at a resonance have small momentum and energy 

relative to the chiral scale. At 12 GeV, the most energetic pions which are produced can 

probably not be treated as soft in a chiral expansion; however, we will see in section 8 

that even at this energy the pion can be treated as soft over a substantial portion of the 

phase space. A naive estimate of the relative rate of BB ± n Zf to BB production would 

say that the first is suppressed by a factor of (p n / An f^) 2 relative to the latter, so that it 

will only be produced significantly at center of mass energy approximately 1 GeV above 

threshold, which is why we investigate not only the resonance regime, but higher center 

of mass energy as well. 
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In our calculation, we need both the chiral and heavy quark approximations to be 
valid. Some care must be taken in the application of the heavy quark effective theory 
to a process with a B and a B meson in the final state. The problem is that the Isgur- 
Wise function (or its analytic continuation) is not useful in the resonance regime, where 
the matrix element (B(p)B(p')\b^ tJ, b\0) varies rapidly as a function of p ■ p' . Because of 
this rapid variation, the form factor is not well described in the resonance region as a 
function of v ■ v', where v and v' are the heavy quark velocities. Furthermore, it is not 
normalized, and drops rapidly to zero beyond the resonance region. This behavior was 
studied by Jaffe in ref. [7], and is expected on the basis of general QCD considerations, 
since the states which can decay into B mesons are not well described as simple Coulomb 
bound states. 

This has several important consequences. First, in the resonance region the calcu- 
lation must be considered as, at best, an order of magnitude estimate. In the language 
of the heavy meson effective theory, this is because there is no well-defined derivative ex- 
pansion; the rapid variation of the cross section with momentum in the resonance region 
is reflected in the heavy meson lagrangian by the presence of higher derivative operators 
"suppressed" only by the QCD scale. This is generally not the case in heavy quark cal- 
culations because of reparameterization invariance[8][9]. However, here, although the 
total cross section for e + e~ — * BBtt respects a reparameterization invariance since the 
cross section depends only on the B meson momenta (p and p') and not separately on 
the heavy quark velocities (v and v') and the residual momenta (k and k'), the rate 
for "decay" of the source is not reparameterization invariant. This is because in the 
effective lagrangian, we work with heavy quark fields at fixed values of v and v', so 
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that we necessarily assume a resonance produced at fixed v ■ v' (not p • p'). Therefore, 
the effective lagrangian must contain nonreparameterization invariant derivative terms 
to reproduce the full cross section, which does respect reparameterization invariance. 
Since the momentum of a B meson is of order \/TbMb, where Tb is the kinetic energy 
of the B meson, the higher derivative terms are always large, even in the heavy quark 
limit so that the derivative expansion is not reliable. 

It is clear that the relative rate predictions for decay to two pseudoscalars, vector 
and pseudoscalar, and two vectors, considered in refs. [10], [11] and [12] are only ap- 
propriate beyond the resonance regime, as has been emphasized by these authors. The 
measured [13] branching ratios for the decay of the -i/>(3S) into D mesons of different 
spins strongly disagree with those predicted by a naive application of the results of the 
heavy quark theory to the resonance regime [14]. The explanation of ref. [15] is that 
the nodes in the ?/>(3S) momentum-space wavefunction result in almost no overlap with 
the final state D mesons except when both are vector particles. In our approach, we 
attribute the large discrepancy between the prediction and the experimental results to 
the presence of higher dimension operators, suppressed only by the QCD scale, which 
give different contributions for decays to heavy mesons of different spins because of the 
significant mass splitting. 

We conclude that the only regime where one would trust the calculation of B 
meson production to better than an order of magnitude is beyond the resonance regime. 
Fortunately, this is the region where the result is most interesting since the rate in the 
resonance regime is too small. When we calculate in the resonance region where heavy 
quark relations are untrustworthy, we view the calculation as a model which should 
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reproduce important qualitative features of the true rate. These include the relevant 
kinematic features, phase space, propagator enhancements, and derivative couplings. 

4. Effective Theory for e+e - BB^* 

We proceed to the construction of the effective lagrangian for BB and BBtt produc- 
tion, keeping in mind the above caveats regarding the application of the heavy meson 
and chiral effective theory to these processes. 

It is useful to divide the range of possible center of mass energies into three regions: 
the resonance regime, above the resonance regime but with the B mesons nonrelativis- 
tic, and high energy, with the B meson fully relativistic. The heavy quark theory is 
straightforward to construct in the third regime, and has been treated in previous work 
[10] [11]. One couples the current to a heavy quark and antiquark of velocities v and 
v'. The heavy quark operator is then evaluated between heavy meson states and the 
appropriate spin symmetry relations between amplitudes can be deduced. 

However, the regions of greatest physical relevance are those at low center of mass 
energy, since at high energy the rate for exclusive production of BB or BBtt is very 
low. Moreover, the experiments of interest are conducted at or near threshold. Since 
the kinetic energies of the B and B are on the order of the QCD scale, both the B 
and B mesons have essentially timelike four velocities in the center of mass frame: 
v» « v't* « (1,0,0,0). 

In the second region — beyond the resonances, but with the B mesons still nonrel- 
ativistic — the effective theory is constructed as at high energy, but with the velocities 
v and v' equal to (1, 0, 0, 0). As before, the virtual photon produces a b and a b quark, 
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each of velocity = (1, 0, 0, 0). In the heavy quark theory, one can couple a source 
to the 6-quark current as b v ' S^j^by (i.e., with the heavy quark spins coupled to the spin 
of the source). When we match onto the heavy meson theory, we can couple the source 
to the mesons as B(f / )S' M 7 At B(z;). This gives the same matrix elements as if we had 
coupled a source to heavy quarks, and then evaluated the heavy quark matrix elements. 

Finally, we consider the resonance regime. If only hard gluons were relevant to 
the binding potential, it would be clear how to construct such an effective theory. The 
matching would again proceed in two steps. First, one would match onto the heavy 
quark theory, and then match the b quark operator onto the heavy meson effective 
lagrangian. However, for a resonance which can decay into B mesons, the binding is 
sufficiently weak that both hard and soft gluons play a role. Hence, it might instead 
be appropriate to match directly onto the low energy heavy quark chiral lagrangian. 
While we still expect that interactions with the light quarks cannot flip a heavy quark 
spin, the interactions between the heavy quarks can flip their spins. The total spin of 
the heavy quarks is conserved, however, and hence the spin of the source is transferred 
entirely to the heavy quark spins. This can be incorporated in the heavy meson theory by 
again coupling the source to the mesons as R(v')S fJ "'y^T5(v) Fortunately, the same heavy 
meson lagrangian describes the matrix elements in both scenarios, since the diagonal 
subgroup of the heavy quark spin symmetry, where both heavy quark spins are rotated 
simultaneously, is sufficient to determine the form of the coupling of the heavy mesons 
to the source. 
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The lagrangian applicable to low-energy production of B and B meson is 
C eS = - ztr{B>K<9 M B a (<;)} - ztr{B>K«9 M B>)} 

+ gtr{B a (v)B b (v)A» bal „ l5 } + gtv{B a (v)B b (v)A^ To } 

+ C S , (4.1) 

£s ~Snr{^B a (v)D: blu B b (v)} 

+ \g'S»tr{^B a (v)A» abl „ l5 B b (v)}. 
Here D = d + V is the chiral covariant derivative incorporating the pion fields. As 

usual, a factor of \[Mb has been absorbed into the heavy meson fields along with the 

position-dependent phase corresponding to the momentum of the heavy quark (so that 

a derivative acting on these fields only gives a factor of the residual momentum), in 

order to suppress the appearance of the heavy quark mass and emphasize the heavy 

quark symmetry. Because this is the low energy theory, no large momentum transfers 

are permitted. At higher energies, the appropriate lagrangian would be the heavy 

meson lagrangian with velocities v ^ v' . The result for two meson production matches 

smoothly, as the difference in residual momenta in the amplitude gets replaced by the 

difference in heavy meson velocities. 

The kinetic and axial coupling terms for the B mesons have been discussed previ- 
ously and result from the straightforward application of heavy quark and chiral effective 
field theories [5]. £5 is the new term and follows from the assumptions described above. 
Note that with the trace the heavy quark spin labels are coupled to 5^7^. 

The coupling A corresponds for a crossed process, e.g., B* — > S7, to the Isgur-Wise 
form factor. Here, we treat the ratio a(e + e~ — > BBjr)/a(e + e~ — > BB) as independent 
of A, although this is strictly true only when beyond the resonance region. (In a full 
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treatment of the resonance region, A would have to be interpreted as a form factor with 
a complicated dependence on the residual momenta, which would not drop out of this 
ratio of cross sections.) 

5. Calculation of a(e+e~ ->• BB^iC*) I o{e+ e~ BB) 

From the lagrangian (4.1) we see that two types of diagrams contribute to BBtc 
production. The pion can be produced "indirectly" by being emitted from a virtual B 
meson through the heavy-meson axial coupling. Or, the pion can be produced "directly" , 
together with the B mesons at a single vertex. This diagram comes from the contact 
term in the lagrangian in which the source couples directly to the B meson and axial 
fields. Examples of both types of diagrams are shown in Figure la. We will see that 
the direct contribution is much the smaller of the two, and so most of our discussion 
concentrates on the indirect contribution. 

In order to compare 

BB ± %t w jth BB as sources of neutral B mesons we normalize 
the BB ± Ti Zf cross sections by dividing them by the cross section for e + e~ — ► neutral B 
mesons, ctq, 



ctq = a(e + e 



B°B° 7 B*°B°, or B*°B*°). 



(5.1) 



We therefore consider the ratios 




-Rpv 



7 ° ± 



BB 



^tt*) + a(e + e 



S*S ± 7T^)}, 



(5.2) 
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and also their sum, 

R = -Rpp + -Rpv + Ryv- (5-3) 

The subscripts denote the heavy-meson content of the BB ± n Zf final state, with a P 
for each pseudoscalar and a V for each vector. In these ratios of cross sections, R a , 
kinematic factors associated with the initial state and also the unknown coupling of the 
source to the B mesons, A, cancel out. 

The cross sections for two-body BB final states are, apart from a common factor, 
each given by the product of a p-wave phase-space factor and a spin-counting factor 
[12]. The total cross section for producing neutral B mesons, 00, is then proportional 
to a sum of such products: 

ao oc rp/p 2 + 4rp/ v 2 + 7^ = P(s). (5.4) 

Here r is the center-of-mass energy that remains after supplying the rest-mass energies 
of the heavy- mesons. For a two-body BB final state, r 3 / 2 oc |AA;| 3 , where Ak is the 
relative three-momentum of the heavy mesons, which is the familiar p-wave phase-space 
factor. The different values of r that correspond to the various BB final states are then 
given by 

rpp =Vs - 2M B 

rpv =Vs - M B - M B * (5.5) 
rw =Vs - 2M B * 

For BB final states, the r Q coincide with the Q values for the decay of the source into 
heavy mesons; however, we will also use the r Q , as defined above, when we consider 
BB ± 7t Zf final states. Note that we are including violation of the heavy-quark symmetry 
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as it enters through the B-B* mass splitting, AM = 46.0 ± 0.6 MeV; this splitting 
significantly affects the available phase space and the off-shellness of intermediate B 



mesons. 



The density of states for a final state of B mesons and a pion simplifies for nonrel- 
ativistic B mesons to 

°° = 256^ ****kS( T -^-^(lAk^?)). (5.6) 



Here kt" = (k + k') fl /2 is the average of the heavy-meson residual momenta and A/c M = 
(k — k'Y is their difference. We have included a factor of (^Mb) 2 for each of the two 
heavy mesons in the final state in order to compensate for the rescaling of the heavy 
meson fields in the amplitudes. If we perform the trivial angular integrals as well as the 
integral over |Afc| (using the 5-function), we obtain the (integrated) density of states, 

D = 1 ^ Mb \Ak\ \k\ 2 d\k\d( cos 9). (5.7) 

Here 9 is the angle between A; and Ak. Then, expressing \k\ and |Afe| in terms of the 
pion energy and momentum, 

p 1T = -2k =>• \k\ 



Ptt_ 

2 



\Ak\= x l4M B (r-E 7V -^ r ) ^2 v /(r-E 7T )M B , 



(5.8) 



the density of states becomes 



D =-^M B \Ak\ \k\ dE 7r d(cos9) 
64tv 6 



(5.9) 



As discussed in the last section, we treat the B mesons as nonrelativistic in the labo- 
ratory frame, working to lowest nonvanishing order in the heavy-meson three-momenta. 

13 



Implicit in our approximations is the realization that the kinetic energy is fairly evenly 
shared between the pion and the heavy mesons: Consider the density of states just above 
(eq. (5.9)). It can be reexpressed in terms of the total kinetic energy, T (= r — m^), 
and the kinetic energy of the pion, T n , as 

D oc T^ 2 (T - T 7r ) 1 ' /2 (T 7r + m,) 1 ^ 

(5.10) 

= T 5/2 x l/2 (1 _ x y/2( x + m J T )l/2 dx 

where x = T^/T is the fraction of the kinetic energy given to the pion. We find that 
as T runs from to ^> m n , the phase-space average of x runs from 1/2 to 4/7, i.e., on 
average the pion kinetic energy about equals the sum of the B meson kinetic energies. 
In the differential cross section, where the density of states is multiplied by the squared 
amplitude, which includes such factors as |Afc| 2 and \p n \ 2 , the powers of x or (1 — x) are 
increased, shifting the average value of x up or down, but the kinetic energy remains 
fairly evenly distributed between the pion and the heavy mesons. 

Therefore, when we work at energies where the heavy mesons are nonrelativistic, 
there is a hierarchy of energy scales, 

T B , E„ (oc M%) < \k B \ (oc M B ' 2 ) < M B , (5.11) 

where Tb is the kinetic energy of the B mesons. Since the dimensionful quantities that 
compensate the different powers of Mb here are of order r, we are essentially working 
to lowest order in r /Mb- Taking sums and differences of the B meson momenta and 
using eq. (5.8), this hierarchy can be reexpressed as 

k°,Ak°, \k\ (oc M%) < |AJfe| (oc M l J 2 ) < Mb- (5.12) 
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In calculating a given amplitude, we retain only the leading term according to this 
hierarchy. Specifically, we drop k and Ak° compared to Mb, incurring errors of order 
r/M B . We also drop \k\ compared to |Afc| and |Afc| compared to M B . This would seem 
to mean dropping terms of order \JVJM~b. However, once an amplitude is squared, 
averaged/summed over initial/final polarizations, and integrated over cos 9, the result 
depends only on fc 2 and Afc 2 (k ■ Ak vanishes in the angular integration). Hence all the 
dropped terms are smaller by a factor of r/M B . 

Later, in section 8, we will consider the contribution to BB ± n^ production from 
a restricted region of phase space where the pion energy is less than a given bound. By 
imposing this cutoff the average value of the pion momentum, and therefore of \k\, will 
be reduced. This only improves our approximation of neglecting \k\ compared to |Afc|. 

We will illustrate these approximations in the simplest case, the calculation of the 
indirect contribution to -Rpp. The two graphs that contribute are shown in Figure lb. 
The graph on the left in Figure lb corresponds to the process S — > BB* — > B(Bir). 
Consider the propagator for the intermediate B* state, which is described by the velocity 
= (1,0,0,0) and a residual momentum = k^ + p% and has a mass that is AM 
greater than that of the B in the final state. The propagator is then 



2(vq + q 2 /2M B - AM) 2(v p n + k ■ p^/M B + m^/2M B - AM) 

*2(v P7T -AM) (5 - 13) 

% 



2^ -AM) 

The first equality comes from the on-shell condition for the B meson, (v + k/M B ) 2 = 1, 
which is just the nonrelativistic formula T B m k 2 /2M B . The next line results from 
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neglecting terms of order k-p^/Ms and m^/M B compared to E n: in accordance with the 
hierarchy in eq. (5.12). The propagator we have obtained, which is inversely proportional 
to the difference in energy between the intermediate and final states, is simply that 
prescribed by nonrelativistic, time-ordered perturbation theory. 

Now consider the part of the amplitude coming from the S — > BB* vertex. From 
our effective lagrangian (4.1) we see that it is proportional to q — k' = Ak +p n - Then 
according to the hierarchy (5.12), we can drop p n = —2k. The amplitude for S — > 
BB* — > B(Btv) is then given by 

Ai(e;v,k,kT) « - 2{E ^f M)u ^v a A^pZe s 1 (5.14) 

where e is the (purely spatial) polarization of the source. The second graph, shown on 
the right in Figure lb, is related to the first by charge conjugation and isospin, so that 
their sum is given by 

A = Ai(e; v, k, k') - Ai(-e; v, k', k) w 2Ai(e; v, k, k') 

(Ak x Pir ) ■ e. 



(En - AM)/, 

Squaring and summing over the averaging over initial polarizations we obtain 

If we substitute for |Afe|, integrate over the density of states, and divide by a"o (not 
including the common kinematical factors corresponding to the initial state) we obtain 
-Rpp. 
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The cross-section ratios for the various BB ± ir Zf final states are given by 

(e»-am) 2 ' a=PP; 

7/4 , 1/2 

(E.-AM) 2 "r" E2_( AM )2 

+ i _i 3/4 a-pv- 

^ (E.+AM) 2 ' " — rv ' 

, + (E T +AM) 2 ' a=VV. 



(5.17) 

Here the z signifies that these are the indirect contributions to R a . The upper integration 
limit is the value of r appropriate to the heavy-meson content of the BB ± ir T final state, 
as given by eq. (5.5). 

At a fixed center of mass energy, the pion attains its maximum energy when the 
heavy meson pair is produced with zero relative momentum (Ak = 0), back to back 
with the pion. From eq. (5.8) we see that, to the accuracy we are are working to, this 
occurs when E n = r a . Hence the upper limit of the phase-space integral over pion 
energy is r a . 

As already discussed, there is also a direct vertex that describes the production of 
B mesons and a pion at a single vertex. For BB ± ir T final states that include at least 
one heavy vector, it is possible to produce the heavy mesons in an s wave. Then, the 
amplitudes are proportional to E^j (Mb fir) while the amplitudes for the indirect graphs 
considered above are proportional to | Afc|p 7r /(Ms/ 7r E' 7r ). At large r, where the cross 
section is dominated by contributions with E n ^> m^, the ratio of direct and indirect 
amplitudes is proportional to E 7T /\ks\ which is of order ^t/Mb- Therefore, the s-wave 
contributions of the direct graphs are suppressed by a factor of r/M B , which is small 
for r inside the region of validity of the chiral expansion. The s-wave direct graphs, 
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Table 2. The ratios a(e + e — > BB ± n Zf )/a(e + e — > £?£?) expressed as percentages. 



v^(MeV) 


n P p 


-Rpy 




R i 


T>d 
"PV 


T>d 

n vv 




R 


10865 


0.12 


0.12 


0.02 


0.26 


0.03 


0.01 


0.04 


0.30 


11020 


0.53 


0.99 


0.70 


2.2 


0.19 


0.10 


0.29 


2.5 


11200 


1.3 


3.1 


3.1 


7.6 


0.59 


0.40 


0.99 


8.6 


11500 


3.3 


9.2 


11. 


24. 


2.0 


1.6 


3.7 


27. 


12000 


8.2 


26. 


36. 


70. 


7.8 


6.7 


14. 


85. 



which cannot interfere with the p-wave indirect graphs, contribute to the R a as 
3^ /2 1 



<{s) = 



^ 2 f 2 P(S) Jm„ 



0, a=PP; 
dE n E 2 P7V (r a - E^) 1 / 2 x { 1/M B , a=PV; 

I /M B , a=VV. 



(5.18) 



The d stands for direct. Note that the source cannot decay into two pseudoscalars in an 
s-wave and a pion without violating either angular momentum or parity conservation. 
Hence Rp P is zero to this order in 1/M B . 

The axial coupling of the heavy mesons, g, has been bounded above by g 2 < 0.5 
[16] using the experimental upper limit for the D* + width [17] which is dominated by 
D* + — > D°n + and D* + — > D + Tt°. Using the maximum allowed value of g 2 and taking 
g' 2 = 1, the numerical values of the ratios R a (expressed as percentages) are given in 
Table 2 for various values of the center-of-mass energy, y/s, and are plotted in Figure 2. 
The sum of the direct contributions to R is also displayed in Figure 2. Because they 
are so small, we neglect them in the rest of our discussion. The overall sum is uncertain 
due to the uncertainties in both the direct and indirect contributions. 

At the T(6S) resonance, the B mesons are produced with a charged pion only a few 
percent as often as they are produced alone. Below this resonance, for example at the 
T(5S), this ratio is negligible, but above it grows with increasing energy as the limited 

18 



three-body phase space for pion production is overcome. For center of mass energy of 
12 GeV, the rate is almost comparable to the rate without a pion. 

Returning to the results for the R l a given in eq. (5.17), we see that in the limit of 
exact heavy quark symmetry, AM = 0, the ratios R a are proportional to one another 

as 

i? PP : R\> Y : R\j Y :: 1 : 4 : 7 (AM = 0). (5.19) 

This limit is approximately realized for AM Cr<C Mb, where the effects of AM ^ 
are minimized and the nonrelativistic treatment still applies. For such large values of r, 
it is also a good approximation to take = 0. Then we find R l a oc r 7 / 2 . The two extra 
powers of r relative to the r 3 / 2 scaling of the BB cross sections reflect the derivative 
coupling and additional phase space factor for the pion. 

The proportions 1:4:7 are the same as was found for the production rates of 
heavy mesons without an accompanying pion in the various spin states (PP, PV and 
VV) using simple spin counting [12] (see eq. (5.4)). We can use the same method to 
understand the persistence of the proportions 1:4:7 when a pion is emitted from one 
of the heavy mesons. 

In S — » BB (pseudoscalars and vectors), the heavy quark spins are fixed so as 
to carry the spin of the source. Charge-conjugation invariance then determines the 
spin state of the light quarks: Consider the BB final state from the point of view of 
heavy quark symmetry, with each heavy meson described by a bispinor. Under charge 
conjugation, the positions and the heavy- and light-quark spins labels of the B and B 
mesons are interchanged. Since the heavy mesons are in a p wave and the heavy quark 
spin state is symmetric (spin one), the light quark spin state must also be symmetric 
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if the final state is to be charge-conjugation odd like the source. Hence the total spin 
of the light quarks must be one. Because the spin of the source is carried by the heavy 
quarks, after the angular positions of the heavy mesons have been integrated over (but 
with the source spin fixed), the total light quark spin points with equal probability in 
all directions. Adding the light and heavy quark spins in each heavy meson to find the 
meson spins, one finds the ratios 1:4:7. 

Now consider the case where a pion is emitted from a heavy meson. Because of 
isospin invariance, we need only consider the case where the pion is neutral and hence 
self-conjugate. Since the 7r° is charge-conjugation even, the total spin of the light quarks 
in the heavy mesons must still be one. After integrating over the angular positions of 
all the particles, the light quark spin distribution will again be isotropic. Hence in the 
heavy quark limit, the heavy meson spin proportionalities do not change when a pion 
is emitted from one of the heavy mesons. 

6. Implications for CP Violation Studies 

We have calculated the BBw branching ratio primarily in order to assess the 
prospects for studying CP violation in the B°-B° system using Yamamoto's pion-tagging 
method [1]. Let us now see if the BB ± n Zf branching ratios obtained in the previous 
section are large enough to produce the statistics required to resolve the small, CP- 
violating asymmetries in the decays of neutral B mesons. The very small values found 
for R at the T(5S) and T(6S) resonances indicate that the pion-tagging method will 
not be useful in the resonance regime, but the rapid increase in the rate with center of 
mass energy suggests that the method could prove useful at higher energies, above the 
resonances. 
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The analysis of this section is due to Yamamoto [1], but uses our results for the 
SS ± 7r =F branching ratios. We follow him in introducing a figure of merit which measures 
the statistical power of a given method of measuring CP violation, 

/ = ae tag d 2 . (6.1) 

It is proportional to the cross section, a, for the process in which the B mesons are 
produced. There is a factor of the tagging efficiency, which is taken to include the 
number of neutral B mesons produced in the process. Finally, it includes the square of 
the "dilution factor" , d, which is given by 

A CP = dsin2/3, (6.2) 

where (3 is a CP-violating angle that appears in the CP- violating asymmetry, A CP , 
measured in neutral B decay. For example, using the pion-tagging method one measures 
the asymmetry 

BB±irT jV(I? — > ipKs) + N(B — > ipKs) l+x 2 M ' v ; 

where N is the number of events from a given decay process, x = 5M/T 0.73 (SM is 
the mass difference between the CP-even and -odd linear combinations of B° and B° 
and T is their common lifetime), and ipKg is a representative CP eigenstate. The utility 
of the figure of merit is that it is inversely proportional to the integrated luminosity 
required in order to measure a CP-violating angle; the larger the figure of merit, the 
easier it is to measure the angle to a given precision. 

In Table 3 we compare the figure of merit for the pion-tagging method, at the 
T(6S) and at 12 GeV, with the familiar methods of generating neutral B mesons for 
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CP-violation studies: T(4S) -> B°B° at an asymmetric collider and T(5S) -> S° J B°7 
at a symmetric collider. Although the asymmetric collider provides the best figure of 
merit, it is worthwhile to determine alternative methods at a symmetric collider in case 
either such a machine is not constructed or the required luminosity is not obtained. The 
relevant benchmark for us is therefore comparison with the figure of merit for studies 
using T — * -B-B7, which are feasible at a symmetric collider. 

The BB(X) cross section is approximately 0.3 nb at the T(5S) and T(6S) reso- 
nances. Using the value of R = 2% at the T(6S) found above, not surprisingly we find 
that the pion-tagging method is not competitive at the T(6S). 

To estimate the BB ± iv Zf cross section at 12 GeV, beyond the resonances, we take 
the (BB)X cross section as 1/2 of its on- resonance value at the T(6S). We further as- 
sume that (BB)X is dominated by pair production of strange and nonstrange B mesons 
and B mesons accompanied by a pion (charged or neutral). We expect that because of 
the sharing of kinetic energy discussed earlier, multipion production will be much smaller 
than single pion production at this energy. To estimate the contribution of BBK and 
BBr] production to (BB)X, we repeated our calculation of BB ± n Zfl production above 
but with K and rj mesons radiated from the heavy mesons, and including SU(3) break- 
ing via the K and r\ masses and decay constants and the B s -B mass splitting. We found 
that BBK and BBr] together are about 1/4 BBn, so we neglect their contribution to 
(BB)X. Estimating B s meson production, as discussed in the following section, by using 
eq. (5.4) but with the values of r a determined by the B s meson masses, we find that 
B s mesons are produced about 80% as often as B°B°. We therefore estimate that the 
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BB ± nT cross section is « 0.15 nb x R/(2.8 + 3R/2) « 0.027 nb. We find a figure of 
merit of 0.005. 

If we extrapolate our results for BB ± TT Zf production to even higher energies, we find 
that at about 12.5 GeV, R is about twice its value at 12 GeV, so that the figure of merit 
is also roughly doubled. We conclude that the figure of merit can be comparable to the 
figure of merit for T(5S) — > B°B°'j at a symmetric collider. Because our calculation 
is not exact, it is conceivable that the figure of merit is even higher. This can only be 
determined experimentally. 

We emphasize that it is unlikely that we have vastly overestimated the rate, since 
it is clear that as one approaches higher center of mass energy that pion production will 
be less suppressed. So long as the center of mass energy is not too high, single pion 
production will dominate. Moreover, we show in section 8 that even if we restrict our 
integral over the phase space to pion energies for which the chiral results are certainly 
reliable, we still get a substantial fraction the total cross section determined from the 
full phase space integration. 

We conclude that Yamamoto's pion-tagging method should be competitive. An 
advantage of running at a higher center of mass energy would be that it could be 
possible to use both BB* production and BB ± n Zf production simultaneously. This 
might augment statistics, and provide a useful check on both methods. 

7. Implications for B s Identification 

Two methods have been used for identifying B s mesons at a resonance: One is 
to scan about the resonance and look for an increased production of strange mesons. 
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Table 3. Comparison of the pion-tagging method with standard 
methods of studying CP violation in the B system. 



Mode 


a 


Stag 


d 


as tag d 2 


B°B° at T(4S) 


0.5 nb 


2 x 0.4 


x/(l + x 2 ) 


0.092 


S°S°7 at T(5S) 


0.05 nb 


2 x 0.4 


2x/(l + x 2 ) 2 


0.015 


BB ± itT at T(6S) 


0.004 nb 


0.8 


x/(l + x 2 ) 


0.0006 


BB ± ttT at 12 GeV 


0.0027 nb 


0.8 


x/(l + x 2 ) 


0.005 



Another method is to identify a B s meson by the lower endpoint of the spectrum of a 
lepton produced in its decay; B s mesons are heavier than nonstrange B mesons and are 
therefore produced with smaller momenta; consequently, a lepton from B s decay will 
have smaller momentum than a lepton from B decay. However, nonstrange B mesons 
produced with a pion also have less momentum than those produced alone, and therefore 
BB-rr is a potential source of background for B s identification. 

We normalize B s production to nonstrange B meson production, 

R s = a(e+e- - B S B S , B S B* S , B*B S , or B*B* s )/a , (7.1) 

where do is the cross section for production of neutral, nonstrange B mesons of spin 
and 1 defined above in eq. (5.1). If we assume that the coupling of the source to the 
heavy mesons is SU(3) symmetric, then up to a common factor both the strange and 
nonstrange cross sections are given by eq. (5.4) but with different values of r a because 
of the B s -B mass difference. (The hyperfine splittings of the strange and nonstrange B 
mesons are essentially equal, however: M B * — M Bs M b * — M B .) To determine the 
importance of the BBtt background we compare R s with 3R/2; R must be multiplied 
by 3/2 in order to include BBtc production. 
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Table 4. Comparison of R s (in %), for various values of Mb s — Mb, with |i? 





80MeV 


105 MeV 


130 MeV 




§#(%) 


T(5S) 


20. 


5.7 


.76 




.39 


T(6S) 


46. 


32. 


20. 




3.3 



Because of the large uncertainty in the experimental value for the strange- 
nonstrange mass splitting, between 80 and 130 MeV [18], we calculate R s for three 
different values of the splitting that span the allowed range: 80, 105, and 130 MeV. The 
results, expressed as percentages, are shown for center of mass energies corresponding 
to the T(5S) and T(6S) resonances in Table 4. 

At the T(5S) resonance, just above B S B S threshold, we find that the B S B S cross 
section depends strongly on the B s -B mass splitting which determines the size of the 
p-wave, B s -B s phase space. As a result, the relative size of the BBtt background also 
depends strongly on this splitting. For the largest allowed value of the B s -B mass 
splitting, the BBtt cross section is about 1/2 of the B S B S cross section while at low 
end it is only a few percent. From the application of heavy-quark symmetry at leading 
order, we expect the B s -B mass difference to be close to the D s -D mass difference. This 
is measured as 99.5 ± 0.6 MeV [6], not far from the center of the allowed range of the 
B s -B splitting. At this point BBtt is about 10% of B S B S . At the T(6S), the B s -B 
mass splitting is less important, and the BBtt background is of order 10% of the B S B S 
cross section. We conclude that the BBtt background is probably not a problem for B s 
identification at either of the T(5S) and T(6S) resonances. 

Of course, it should be borne in mind that this calculation assumed SU(3) symmetry 
and was based on heavy meson effective chiral theory applied in the resonance regime. 
However, this naive calculation indicates that BBtt background should not be a problem. 

25 



8. Measurement of g 

So far, we have integrated over the entire phase space because both CP violation 
studies and the identification of B s mesons rely on the total rate. In this section, we 
instead focus on the region of phase space where we expect the calculation to be reliable. 
This serves two purposes. First, it allows for the possible extraction of the axial coupling 
constant of the B mesons, g. It would be useful to directly extract g in this way and 
compare to the bounds on g in the D system. This method is probably not useful for 
extracting g in the D system however as we show below. Second, we can establish a 
reliable lower limit on the branching fraction to BBtt for center of mass energies beyond 
the resonance regime. 

Because we want a reliable prediction, we focus on energies beyond the resonance 
region. We also need to determine when the derivative expansion of the heavy quark 
chiral lagrangian is sufficiently reliable that we can trust the leading order (in deriva- 
tives) result. As stated in ref. [5], we need both v-p^ and v' -p n to be small. How small 
depends on the cutoff for the theory; this may be the chiral symmetry breaking scale, 
A x or it may be smaller (see ref. [19] for a discussion). We call the cutoff A. 

Notice that in order for both v-p n and v' -p n to be small, v-v' and v° are restricted. 
We see this by adding both constraints together, which yields iPE^ < A. The lowest 
possible E n is m n , implying v° < A/m f which in turn implies v ■ v' < 2(A/m 7r ) 2 — 1. 

This means that even for the pion emitted at threshold, the momentum of the 
heavy mesons is constrained if the chiral expansion is to be valid, and furthermore, 
that the best region to apply the heavy meson lagrangian will be in the region where 
U° 1 but above the resonance region, where the calculation will be valid over the 
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maximum possible range of pion energies. Therefore, we will concentrate on extracting 
g in the regime where the B mesons are nonrelativistic. This is clearly the best place 
from an experimental vantage point, and as we have argued, is probably also the region 
where the branching fraction to B mesons and a single pion is maximal. One might 
also hope to extract g for D mesons. However, at CLEO, where D mesons are copiously 
produced, v° D 3. Therefore, even with a high cutoff for A of order 1 GeV, one could 
only integrate to pion energies less than 2m^. Even with this limited region of phase 
space, v ■ and v' ■ p^ are very close to the cutoff so that the extraction of g for D 
mesons is probably not reliable. So we focus on the extraction of g in the B system. 

Recall that our expression for R — the ratio of the cross section for production 
of the various BB ± rc T final states divided by the cross section for BB production - 
was expressed as an integral over pion energies (see eq. (5.17)). We now integrate the 
differential cross section over only up to some maximum pion energy, i?™ ax . We will 
refer to the result as R(s; E™ &x ). Ultimately one wants to choose _E™ ax as the maximum 
pion energy for which we expect the cross section to be reliable. The allowed range of 
^imax j g De t wee n and rpp. 

In Figure 3, we plot R(s; £™ ax ) as a function of £™ ax for yfi of 11.5 GeV, 11.75 GeV 
and 12 GeV. Note that for £™ ax < 500 MeV, R(s; £™ ax ) is roughly independent of y/s 
in this range. This is because at these energies, the p-wave phase-space factor for the 
heavy mesons is not much different for heavy mesons produced with a soft pion than for 
heavy mesons produced alone, and approximately cancels out in the ratio i?(s;£'™ ax ). 
Accordingly, for £™ ax = 400 MeV we find that R(s; £™ ax ) is roughly 10% for any value 
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of > 11.5 GeV; for £™ ax = 500 MeV we find R(s; -E™ ax ) slowly varies from about 
15% to 20% as v/s runs from 11.5 to 12 GeV. 

For larger values of i?™ ax , the dependence of R(s; i?™ ax ) cannot be ignored. If we 
allow a large value of the cutoff, £;™ ax = 800 GeV, we find that we gain rather little 
at v^i = 11.5 GeV, moving up to a value of R(s; £™ ax ) just above 20%, whereas at 
11.75 GeV we have increased to over 35%, and at 12 GeV to over 45%. 

This analysis also allows us to conclude that we can reliably predict a large value of 
_R a , since within the regime of reliability of our calculation, the ratio R a is very large, 
greater than 45% at y/s = 12 GeV if we allow a cutoff of 800 MeV. It is encouraging 
that this rate is so large. 

It is important to recognize that one can first test that the heavy quark symmetry 
relations apply, by studying the relations among pseudoscalar and vector B meson 
production without a pion. If these prove valid, one can then proceed to measure the 
cross section with a sufficiently small cut on the pion energy so that the chiral theory 
is applicable. By varying the cutoff, one could test where the results disagree with our 
prediction, indicating that the derivative expansion of the heavy meson effective theory 
has broken down. This should permit a reliable extraction of g. 

9. Conclusion 

It is clear that the proposal of Yamamoto is quite interesting. It appears that there 
might be a sufficiently large rate for self-tagging B meson events for this to be a viable 
method of studying CP violation for neutral mesons at a symmetric collider. Despite the 
limitations of our calculation, we can nevertheless establish several interesting results. 
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Within the resonance regime, the process we consider will probably not occur at a 
sufficiently large rate to compete with the more conventional proposal for the study of 
CP violation at a symmetric collider, namely BB* production followed by B* — > S7. 
At center of mass energy of about 12000 MeV, pion emission from a B meson pair could 
occur as often as not. Although the calculation here is reliable only over about half the 
range of pion energies, it is clear that the rate for a single pion accompanying the B 
mesons is large, even from this restricted phase space. 

It could therefore be useful to run at center of mass energy above the resonance 
region. In this region, one has the advantage that BB* and BB ± ix^ production should 
both be large. With two different methods of looking for CP violation, more reliable 
results might be obtained. And of course, the BB ± 7r Zf process has a much cleaner tag 
and employs known technology. 

Furthermore, it is likely that we have been conservative in our evaluation of the 
utility of Yamamoto's method. We assumed a drop in cross section by a factor of 2 at 
center of mass energy 12 GeV. In addition, we have neglected production of resonant 
B meson states which would decay primarily to Bn. Because the production of such a 
resonance with a low-lying B meson state is not p-wave suppressed, there could be a 
larger rate for BB ±r K^ than we computed. 

From our calculation, we have also established that at the T(5S), the production 
of BBn is probably not a large background. 

Finally, we have shown that the rate for BB ± 7r T outside the resonance regime 
should allow for a reliable extraction of the axial pion heavy meson coupling constant, 
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g. In general, running above resonance but in the regime where B mesons are nonrela- 
tivistic could be the best place from the point of view of testing heavy quark relations 
at energies sufficiently high that they should be reliable, but not so large that exclusive 
modes are suppressed. It is encouraging that the rate is largest in the range of energies 
where the calculation should be most reliable. 
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Figure Captions 



Figure 1. Diagrams that contribute to BBtt production. The blob represents the source S^. 

(a) Examples of the two kinds of graph: indirect and direct, (b) The two indirect 
graphs that contribute to the production of two pseudoscalar B mesons and a pion. 

Figure 2. The ratio, R : of the SS ± 7r =F and B°B° cross sections as a function of the center 
of mass energy. Both indirect and direct contributions are shown. Also shown are 
the individual indirect contributions, R l a , coming from specific BB ± ir Zfl final states. 
The couplings were taken as g 2 = 0.5 and g' 2 = 1. 

Figure 3. The (indirect) contribution to R for pion energies less than a cutoff E™ ax for center 
of mass energies of 11.5, 11.75 and 12 GeV. 
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